Attitude Determination and Prediction 
of Spin- Stabilized Satellites 
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(Manuscript received February 17, 1964) 

Techniques for both altitude determination and prediction for spin- 
stabilized satellites are developed. Their use is demonstrated using Tel star 
I and II satellite data. It is shown that an inclined dipole model of the 
earth's magnetic field and the method of averaging the gravitational and 
magnetic torques over each anomalistic period of the satellite permits attitude 
predictions to within a few tenths of a degree of determined values in most 
instances. In those few cases where departures are above one degree, explana- 
tions are presented to shoiv the reason for such discrepancies. 

The usefulness of combining optical flash and solar sensor data for 
attitude determination and their inherent accuracy are demonstrated. Opti- 
cal flash data can provide loci with a resolution of 0.1°. /Solar sensor loci 
are resolved to within 1°. 

All techniques described have been consolidated into working computer 
programs which follow closely the mathematical analysis presented. A num- 
ber of important supporting calculations such as the solar position, sidereal 
time, orbit updating, etc. are also developed. Because of the complexities of 
the mean torque and gyroscopic equations, the precessional techniques 
presented are most useful in computer embodiments. 
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I. INTRODUCTION 

To maintain a defined attitude in space, the Telstar I and II satellites 
were spin stabilized. By this method of passive attitude control, a 
satellite is rotated about an axis of symmetry and consequently exhibits 
the characteristics of a gyroscope. In the absence of disturbing torques, 
the satellite's spin axis maintains its spatial orientation fixed with 
respect to an inertial reference frame throughout its orbit. For the 
Telstar I and II satellites, this is desirable because of certain required 
attitude constraints. First of all, the satellite communication antenna 
is not omnidirectional. More energy is radiated along the equator of 
the satellite than along its spin poles, as shown by the antenna pattern 
of Fig. 1. This fact dictates an attitude for which the line of sight from 




Fig. 1 — Antenna pattern at 4170 inc. 
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Fig. 2 — Skin temperature distribution vs polar angle. 

a ground station to the satellite avoids its poles. The second attitude 
constraint involves temperature considerations. A good degree of tem- 
perature control is obtained by orienting the spin axis so that it is 
nearly perpendicular to the satellite-sun line. In this manner, tempera- 
ture balance is maintained by the satellite's spin as indicated in Fig. 2. 
Here, <p is the solar offset angle, defined as the angular departure of the 
spin axis from perpendicularity with the satellite-sun line. Solar offsets 
of about 15° result in temperature deviations of about 150° (see Fig. 2) 
and are tolerable. 
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Aii ideal orientation from a communications standpoint would be to 
have the spin axis nearly parallel to the earth's surface as it passes over 
any ground station. A spin axis perpendicular to the orbital plane would 
accomplish this, but would produce a maximum axis tilt toward the sun 
equal to the sum of the orbital inclination and the earth's 23$° tilt with 
respect to the plane of the ecliptic (see Fig. 3). Under these conditions 
the spin axis of the first Telstar satellite would have a maximum solar 
tilt of 08°. This exceeds the 15° tilt limit dictated by temperature 
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Fig. 3 — Geometry of spin axis perpendicular to orbit plane. 

balance. On the other hand, if the spin axis is made perpendicular to the 
ecliptic plane (which automatically insures its perpendicularity to the 
sun), the axis tilt with respect to the line of sight from any ground 
station over which it may be passing will range from 90° to 90° minus 
the sum of the orbit inclination and ecliptic obliquity or 22° (see Fig. 
4).* From a communications standpoint this is tolerable, since nulls 
in the antenna pattern arc major only within about 15° of the spin 

* With this orientation, stations south of the satellite's instantaneous earth 
latitude would experience angles less than 22°, such as station B of Fig. 4. The 
major ground stations for Telstar I and II, however, are all at latitudes above 
43°. Since the inclinations of the Telstar I and II satellites are 45° and 42.7° re- 
spectively, these stations are almost always north of the satellite. 
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Fig. I — Geometry of spin axis perpendicular to ecliptic. 

poles (see Fig. 1). Rising and setting satellites may experience a greater 
tilt to the ground station line of sight, but this is unavoidable. 

Because of these considerations, the Telstar I and II satellites were 
launched with their spin axes as nearly perpendicular to the ecliptic as 
the powered flight trajectory of the Thor-Delta launch vehicle would 
permit and still meet certain orbital requirements, such as inclination, 
apogee height, and perigee height. 1 The predicted attitudes of these 
satellites at orbit injection are calculated from telemetry data from the 
first and second rocket stages. Table I lists these attitudes as the right 
ascension and declination direction of the north pole of the satellite. 



Table I — Attitude of the Telstar I and II Satellites 
at Injection, Determined from Powered Flight Data 





Initial Attitude 




Right Ascension 


Declination 


Telstur I 
TelBtar II 


83.73° 
82.23° 


-60.80° 
-57.31° 
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The north polo is defined here as the direction of advance of a right- 
handed screw turning in the direction of the satellite's spin. It also 
happens to be the spin pole which carries a helical telemetry antenna. 

The initial attitudes of the Telstar I and II satellites are changed 
by gyroscopic precession. This motion of spin-stabilized satellites is 
chiefly produced by both magnetic and gravitational torques. The 
former is a result of interaction between residual and eddy-current- 
produced magnetic fields of the satellite and the earth's magnetic field. 
The latter is produced by differential gravitational forces acting across 
the body of the satellite. In the present cases, the magnetic torques are 
several orders of magnitude in excess of the gravitational torques. 

It is necessary to predict precession to allow proper scheduling and 
planning for satellite use, to resolve certain attitude determination 
ambiguities, and to sensibly plan alteration of satellite attitude in a 
prescribed manner when needed. (See Section V for a description of this 
technique.) It is the purpose of this paper to outline the methods of 
attitude determination and correction, to develop the precessional 
theory, and to show the application of these to the Telstar I and II 
satellites. 

II. ATTITUDE DETERMINATION — GENERAL REMARKS 

The attitude of the Telstar I satellite has been determined through 
the analysis of two sets of data: the time of optical flashes of sunlight 
from three mirrors attached to the surface of Telstar and the current 
produced by six onboard solar sensors located on the ends of three 
orthogonal axes. 

The first of these sets of data in combination with the spatial position 
of the satellite, sun, and observer's position determines a locus of possible 
spin axis positions which would result in the observed flash. This locus 
describes a cone in space whose axis is the mirror normal. The six solar 
sensors are designed to determine the direction of the satellite- sun line 
with respect to a satellite frame of reference and thus to the spin axis. 
Knowing this angle, again there is defined a conical array of possible 
spatial spin axis directions. The intersection of the optical cone with 
this solar sensor cone should, therefore, determine two possible attitudes. 
A priori knowledge of the approximate attitude as provided by launch 
data and/or a succession of measurements and predictions over an 
interval of several days permits the determination of a unique spin 
axis direction in space. 

In practice, the solar sensors present a few difficulties. First of all, 
the deduction of the sun's position from solar cell current entails inferring 
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a solar angle between each cell's normal and the sun's direction. Such a 
direction from each of three cells uniquely establishes this solar aspect.* 
Thus calibration curves relating electrical output to light intensity must 
be employed for each cell along with temperature corrections. While 
D. W. Hill 2 has prepared a computer program to take the drudgery out 
of this work, it remains difficult to calculate the solar angle for a cell if* 
it is illuminated in a direction far from its normal. Moreover, the solar 
sensor data, which are reported every minute by telemetry, were found 
to be not always mutually consistent. Often, the solar direction calcu- 
lated over 30 minutes from a succession of telemetry frames had a spread 
ranging from 1 ° to over 8°. A correlation between these deviations and 
the spatial position of the satellite exists which suggests a biasing of 
solar cell data, on occasion, by the earth's reflectivity. 

To show the geometry of this situation, consider the satellite position 
shown in Fig. 5. Here one half of the satellite to the right of line AB is 
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Fig. 5 — Satellite shadow geometry. 

illuminated by the sun. At least one half, to the upper left of line CD, 
is illuminated by the earth in reflecting about 38 per centf of the sun's 
total incident light. Thus a solar cell in region BOD would record only 
direct sunlight, one in region AOC would record only earthshine, and 

* The solar aspect is the angle between the satellite's spin axis and the satel- 
lite-sun line. 

t This is not t lie albedo but the earth's mean reflectivity or ratio of mean earth 
brightness measured at a spot along the earth-sun line relative to the brightness 
of a perfectly diffusing disk of the same size and at the same distance that the 
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one in region AOD would record both. There would be no light from 
either the sun or the earth in region COB. By the difference in solar 
illumination and earthshine, data from cells in regions BOD and AOC 
are easily separated and the latter disregarded. However, those cells 
in region AOD measure both sunlight and earthshine in an inseparable 
manner and therefore report erroneously the sun's position. If AOD is a 
large angle, the spin of the satellite will carry solar sensors into that 
region frequently and result in sizable variations in successive solar 
aspect determinations. Since AOD equals SET, as the subsatellite 
position (T) departs from the subsolar point (S), the solar cell data 
must be carefully interpreted to avoid false conclusions. Under these 
conditions, it would be wrong, for example, to simply calculate the mean 
of all solar aspect determinations over a succession of telemetry frames. 
This procedure, in general, would not yield a good estimate of the true 
aspect. A better technique would be to calculate the mean for only 
those solar aspect determinations in which all three cells entering into 
the determination exhibited current readings above those which could 
possibly be produced by either the earth's reflectivity or low angles of 
solar illumination. Operating in this fashion, the true solar aspect can 
be determined to within about 1°. 

Since the maximum attitude resolution obtainable using optical 
flashes from the Telstar satellites is about 0.1°, it was decided to rely 
on these for attitude determinations insofar as possible. Two groups 
of flashes close together in time are needed for an attitude fix, however, 
and in cases where only a single one existed, the attitude was determined 
by a combination of mirror flashes and the solar sensor data previously 
described. 

The optical reflections are characterized by a series of intermittent 
flashes provided by the spin of the Telstar satellites. The time midpoint 
of these flash series is determined by photoelectric equipment 3 at Bell 
Telephone Laboratories in Holmdel, N. J. The time of each flash series 
determines a conical locus of possible spin axis positions about the 
flashing mirror's normal. The tip of the spin axis vector, therefore, 
lies on a circle on the celestial sphere. One of the intersections of two 
such circles defines the attitude, provided the two corresponding flash 
series occur close enough together in time so that no appreciable preces- 
sion occurs during the separation interval. Since the mirrors employed 
have their normals far removed from the spin axis (68° and 95°), it is 

earth is from the given spot. Quite naturally a phase law applies which reduces 
the light reflected to points off the earth-sun line. Off-line brightness vanes 
crudely as the ratio of observable illuminated area of the earth as seen from the 
point in question to the total observable illuminated area seen from a point on 
the earth-sun line. 
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always possible to select the valid locus intersection from past per- 
formance and attitude predictions. 

III. ATTITUDE DETERMINATION — COMPUTATIONAL TECHNIQUE 

A computational technique has been developed for the calculation 
of spin loci. It consists of four basic parts, which are: 

(i) the determination of the right ascension and declination of the 
normal of the Hashing mirror; 

(ii) the determination of the right ascension and declination of the 
sun in the case of solar cell data; 

(Hi) the construction of a circular locus on the celestial sphere cen- 
tered on the right ascension and declination of the mirror normal and 
having a radius of (>8° or 95°, depending upon the mirror involved (in 
the case of solar aspect data from the sensors, the locus is centered on the 
right ascension and declination of the sun and has a radius equal to the 
measured solar aspect angle) ; 

(iv) the plotting, by computer microfilm techniques, of these circles 
and others determined from additional flash series and solar aspect data 
to ascertain intersections and corresponding satellite attitudes. 

3.1 The Right Ascension and Declination of the Mirror Normal 

Consider the generalized minor orientation shown in Fig. fi. Here, a 
ray from the sun strikes the mirror and is reflected toward a particular 
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ground station. For this to be possible, the mirror normal must bisect 
the angle 2a and lie in the plane determined by the mirror-sun line and 
the mirror-ground station line. 

Sketched around the mirror in Fig. 6 is the celestial sphere. On this 
sphere the solar direction, as seen from the mirror, is indicated by its 
right ascension, a s , and its declination, 8 a . In like fashion, the direc- 
tion of the satellite mirror as seen from the ground station is specified 
as a,8 and the direction of the mirror normal is a.,8,, . 

To determine a„ and 8„ , we begin by solving for arc g, and A in the 
spherical triangle 1, 2, 3: 

cos g, = sin 8 sin 8 S + cos 8 cos 5 6 cos (a„ — a) ( 1 ) 

. sin 8 S — sin 8 cos g, , -, 

cos A = ^^ \*i 

cos 5 sin g 8 

where angles g s and A may have values from 0° to 180°. 

By triangle 1, 2, 4, one obtains the declination of the normal as 

8„ = sin -1 [sin 5 cos g s + cos 8 sin g„ cos A ] (3) 

where 8„ may have values between ±90°. 
Using the same triangle we may write 



C = sin -1 



sin (g „ + a) sin 8 n l ,^ 

sin A J 



where 

A = 18 ° - °- . (5) 

C may have values between 0° and 180°. Since a„ may either exceed a 
by C or be less by the same amount owing to the two possible orienta- 
tions of triangle 1, 2, 4, we have 

a„ = a + SIGNF (C,a s - a) (6) 

where SIGXF, a common computer symbol, indicates that the algebraic 
sign to be affixed to C shall be determined by the quantity a, — a. 

Right ascension is measured eastward from the vernal equinox as a 
reference through 300°, and if a, and a should lie on opposite sides in 
this reference (C remaining less than 180°), (6) becomes 

a„ = a + SIGNF (C,a - a g ) . (7) 

Thus (6) is employed if | a. — a \ - 180° is negative and (7) if other- 
wise. 
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3.2 The Right Ascension and Declination of the Sun 

The apparent right ascension and declination of the sun at any 
specified time is computed from the mean orbital elements of the sun." 1 
These may be expressed beginning with the true solar mean anomaly in 
degrees as 

.1/, = 808.47583 + 0.9850002070r/' - 0.000150T 2 - 0.000003T 3 

where 

T = the time in Julian centuries of 36525 ephemeris days from 

January 0.5, 1900, ephemeris time 5 
d' = ephemeris days from same epoch. 

We may express the above equation in a more useful form for the present 
calculations by changing the epoch to January 1.0, 1900 and writing 
an equivalent expression as 

M, = 357.41283 + 0.985000207tf 

where 

d = ephemeris days since the 1900 epoch 

= mean solar days since 1900 epoch + 1 second for the years 
1903, 1904, 1965. 

One obtains the apparent mean solar anomaly, used in the present calcu- 
lations to determine the apparent position of the sun, by antedating for 
the solar light transit time. For this reason d is increased by 0.005375 
day, which is the light transit time at mean solar distance. Since the 
earth is about 3 X 10'' miles closer to the sun in winter as compared to 
summer, this can produce an error in apparent solar position of about 2 
seconds of arc. 

The eccentricity of the earth's orbit is 

e = 0.01075104 - 0.000041807 1 - 1.20 X 10~T 2 

- 0.01700254 on January 1.0, 1900. 

Also, the degrees of mean celestial longitude of the perigee of the sun's 
mean orbit about the earth as a reference is 

L = 279.09008 + 0.9850473354^/' + 0.000303 T 2 
= 282.25247 + 0.470084 X 10~ 4 rf. 
Finally, the mean obliquity of the ecliptic in degrees is 

c = 23.452294 - 0.01301257' - 1.04 X 10~ 6 7 t2 + 5.03 X 10~ 7 T 3 . 
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The apparent mean anomaly and eccentricity along with Kepler's 
equation permit the calculations of the apparent true anomaly of the 
sun by standard techniques. '"' This along with the apparent mean longi- 
tude of solar perigee and mean obliquity of the ecliptic allows a direct 
determination of the apparent right ascension and declination of the 
sun at any date by the simple geometry shown in Fig. 7. 

The apparent right ascension and declination of the sun computed 
from its mean elements in this manner can depart from the tabular 
values (which include nutation) of The American Ephemeris and Nauti- 
cal Almanac for the years 1962, 1963, and 1964 by about 5 seconds of 
arc as a maximum. This is well within the approximate 1° error of the 
measuring techniques for solar aspect as well as the 0.1° for the mirror 
technique. 

3.3 Construction of the Locus Circle 

As previously stated, the direction of the spin axis may be specified 
as a circle on the celestial sphere, centered on the right ascension and 
declination of the mirror normal. This locus may be generated by using 
the spherical triangle shown in Fig. 8. Here the dotted curve indicates 
the attitude locus of all possible positions of the spin axis. Let aiSi be 
any point on this circle making a fixed arc F with a„ , 8„ . In the case 
of cither the Telstar I or II satellite, F will equal 68° or 95°, the angles 
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Fig. 7 — Geometry of the sun's position. 
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Fig. 8 — The attitude locus, 
the mirror normals make with the spin axis. From Fig. 8, 



5* = sin [sin 5„ cos F + cos 5„ sin F cos .4] 



(8) 



where 



A = a running variable which takes on values from to 360° to 

generate the locus 
d t = declination of a locus point which ranges from —90° to +90°, 

in general. 



Also 



sin B = 
cos B — 



B = tan 



sin F sin .4 
cos 5; 

(cos F — sin 8 n ) sin 5/ 
cos 6 n cos 5; 

i sin B 



cos B 



and 



on = a„ + B. 



(9) 
(9a) 
(10) 

(10a) 



In general, B may range from 0° to 360°. All quadrant ambiguities 
presented by (10) are resolved by noting the algebraic signs of (9) and 
(9a). 

If there is a time uncertainty in measuring the midpoint of a flash 
series, this will result in more than one possible mirror normal, a„5„ , 
and hence a number of attitude loci, since the satellite and sun will 
occupy successive positions along their paths within the time error. 
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If during a single sky trajectory a ground station records two separate 
flash series (or a single series plus solar aspect data), each will generate 
a circular locus. The intersections of these loci determine attitude with 
an ambiguity of two. That is, there are two intersections and therefore 
two possible satellite attitudes which satisfy the two flash series. A 
priori knowledge of the attitude from previous measurements or initially 
from launch parameters together with some knowledge of the expected 
precession will permit selecting the proper intersection. Knowledge of 
expected precession is also needed to determine whether or not the flash 
series arc close enough together in time to neglect incremental precession 
during the time between the series. Such precession can alter the position 
of the loci intersections. 

Fig. 9 shows an example of the attitude determination of the Telstar 
II satellite for passes 135, 136, 199, 272, and 472. Solid lines indicate 
attitude loci determined from mirror flashes. The line of zeros is typical 
of those determined by solar sensor data. Two attitude determinations 
are shown by intersections on passes 135, 130, and on pass 272. 

IV. ATTITUDE PREDICTION 

Acting upon an orbiting spin-stabilized satellite to produce precession 
are certain disturbing torques. Those to be considered here are gravita- 
tional and magnetic torques. In the case of the Telstar I and II satellites, 
these are dominant over atmospheric drag torques, solar radiation 
torques, electrostatic torques, and others. 

4.1 Assumptions for Gravity Torque Calculations 

Differential gravity forces acting across the body of the satellite can 
produce torques which tend to rotate the body. These forces exist 
simply because the strength of the earth's gravity field is a function of 
the distance from earth. A body of finite size must, therefore, experience 
such torques. 

For the purpose of calculating the mean gravity torques, the earth 
is assumed spherical with its radius equal to the equatorial radius. For 
the Telstar satellites, gravity torques are at least an order of magnitude 
less than magnetic torques, and therefore neglecting the earth's oblate- 
ness produces at worst only a second-order error. 

It is also assumed that the moments of inertia of the satellites about 
all axes perpendicular to the spin axis and passing through the mass 
center are equal. The orbit is assumed elliptical and Keplerian, since 
the earth's figure produces but second-order effects over a single satel- 
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Fig. 'J — Computer plot of mirror and solar loci. 

lite period. The orbit, however, for each succession of calculations will 
he updated using mean orbital elements, which include all secular per- 
turbations, to produce the Keplerian orbit of the best fit at the beginning 
of each calculation. 



4.2 Assumptions for Magnetic Torque Calculations 

In order to construct a mathematical model which on one hand rea- 
sonably well represents the physics of the situation, but on the other 
hand does not, by its complexity, produce numerical equations costly 
to compute, the following simplifying assumptions are made: 
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(i) The magnetic field of the earth is represented by a dipole centered 
at the earth's center and inclined degrees (equal 11.4°) to the earth's 
spin axis. 7 •" 

(ii) Mass symmetry about the satellite's spin axis exists, as in Sec- 
tion 4.1. 

(Hi) The effects of magnetic moments transverse to the spin axis are 
either negligible or average out due to the satellite's spin. 

(iv) For the calculation of mean or net magnetic torques over an 
orbital period, a Keplerian orbit is assumed as in Section 4.1. 

4.3 Assumptions for Gyroscopic Equations of Motion 

It will be useful in simplifying the gyroscopic equations of motion 
(see Section 4.7) to assume that the satellite angular momentum vector 
coincides with the spin axis and equals the spin rate times the moment 
of inertia about the spin axis. This is the same as assuming the entire 
angular momentum of the satellite is the result of its spin alone, and 
neglects that small amount provided by the precession of the spin axis 
itself. Precessional dampers 9 are provided on the Telstar satellites to 
prevent coning or rapid changes in attitude at rates comparable to the 
spin rate. This, more than ever, makes the assumption quite reasonable. 

4.4 Coordinate Systems 

It will be convenient to establish certain useful coordinate systems 
and their interrelationships prior to the torque calculations. These will 
be defined and related by Euler-type axis rotations expressed by matrices. 
All coordinate systems are rectangular and right-handed in the con- 
ventional sense such that the rotation of an X axis into a Y axis de- 
termines the positive direction of a Z axis as the direction of progress of 
a right-handed screw. Each system will be named, described, and inter- 
related in that order. 

4.4.1 The Earth-Centered Inertial System (IS) 

In this system, the three mutually perpendicular X, Y, and Z axes 
have their common origin at the earth's center. X contains the vernal 
equinox and increases positively from the earth's center in its direction. 
Y is perpendicular to X in the earth's equatorial plane. Z contains the 
earth's spin axis and increases positively toward the north celestial 
pole. This is the basic system to which others will be referred. 



ATTITUDE DETERMINATION AND PREDICTION 



1673 



4.4.2 Rigid-Body Systems 

Two rigid-body systems will he used. Both have their origin at the 
satellite mass center. The first of these is to be known as the SANOR 
or satellite nonrotating system. In this system, the z axis defines the 
satellite's spin axis, being positive in the direction of advance of a right- 
hand screw spimiing with the satellite. The x axis is in the satellite 




TO VERNAL 
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Fig. 10 — The IS, SANOR and SAR systems. 

equatorial plane, \p degrees rotated from X (sec Fig. 10) about Z. This 
axis defines the xy-XY plane intersection. The y axis is in the satellite 
equatorial plane orthogonal to x. At this point, we note that the rota- 
tional equations of motion for the satellite will be developed in this 
SANOR system. 

The SAR or satellite rotating system, as its name implies, differs from 
the SANOR system in that it rotates with the satellite. It has axes 
x', y', z', where z' coincides with z, and x' and y' are defined as being 
rotated <p degrees from x and y respectively about z. 

To relate the SANOR and SAR systems to IS, we proceed as follows 
(refer to Fig. 10). The ^ rotation about Z yields the following matrix 
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I) = 



cos \p sill \f/ 

—sin \f/ cos \p 

1 



(11) 



The 9 rotation about x yields 



C = 



The <p rotation about z gives 



B = 



1 

cos 6 sin 6 
— sin 6 cos 6 

cos <p sm<p 

— sin ^> cos ip 

1 



(12) 



(13) 



X 




X 


y 


= (CD) 


Y 


z 




Z 



To transfer points in the IS to points in the SANOR system, we have 

(14) 
Z 

or, just to shorten the notation 

x = (CD)X (15) 

where 

(CD) = multiplication of the C and D matrices. 
In like fashion 

a' = (BCD)X (16) 

= AX (17) 

where 

A = (BCD) = multiplication of the B, C, D matrices. 
Carrying out this multiplication, we may express .4 as the following 
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A = 



cos <p sin \p sin tp sin 6 

+ cos 8 cos \f/ sin <p 
— sin <p sin \p cos ^ sin 8 

+ COS COS \f/ COS y> 

— sin 6 cos ^ cos 



cos ip cos ^ 

— cos 8 sin ^5 sin \p 
— sin v 3 cos \f/ 

— cos sin ^ cos <p 
sin sin i/' 

The matrix CD may be obtained from A by letting <p = 0, 

cos ^ sin i/' 

( CD ) = —sin yp cos 8 cos ^ cos sin 8 
sin i/< sin 8 —cos i/' sin cos 8 



(18) 



(19) 



Quite obviously the inverse operations apply by taking the transposes 
of the matrices. Indicating a transposed matrix by the symbol ~, we 
have 



and 



X = (DC)x 



X = Ax'. 



20 



(21 



4.4.3 Orbital Coordinate Systems 

Let an orbit defining system (ORDEF) be described along the lines 
shown in Fig. 11. Here x„ is the intersection of the orbit plane and the 
earth's equator plane with +.<•„ drawn toward the ascending node of the 
orbit from the earth's center. Axis z„ is normal to the orbit plane, posi- 
tive in a direction a right-handed screw along oz g would advance if 
turned in the direction of the satellite's orbital motion. Axis //„ completes 
a right-handed system by being mutually perpendicular to .r„ and z„ . 

The satellite defining system (SADEF) has its z s axis collinear with 
z„ , but x» passes through the instantaneous satellite position, w degrees 
from x„ . Axis x s therefore will be referred to as the local vertical of the 
satellite. 

The ORDEF and SADEF systems are related to the inertial system 
(IS) in a manner strictly analogous to the SANOR and the SAR mat- 
rices described in the previous section. The only change in the mat- 
rices is that 
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ft replaces \p 
i replaces 6 
o) replaces <p. 
Thus we may write 

x = (CD)iX (22) 

x = {DG)a. (23) 

x, = AiX (24) 

X = AiX. (25) 

where the subscript i indicates the ft, i, w substitution. 

4.4.4 Magnetic Coordinate System 

A magnetic coordinate system (MAG) based on an inclined dipole 
model of the earth's magnetic field is constructed as follows. Let axis 
x m define the intersection of the geomagnetic equatorial plane with the 
earth's geographical equatorial plane, positive toward the ascending 
node of the geomagnetic equator (17 from X, see Fig. 12). Axis z m is 




Fig. 11 — The ORDEF and SADEF coordinate systems. 
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CELESTIAL 
EQUATOR 



GEOMAGNETIC 
EQUATOR 



Fig. 12 — The geomagnetic coordinate system. 

normal to the geomagnetic equator plane, positive toward the north 
geographical hemisphere, and 0° from Z. Here /3 is simply the magnetic 
dipole inclination to the earth's spin axis. Axis y m is in the geomagnetic 
equatorial plane orthogonal to x m . 

Transformations from the MAG system into the IS system proceed 
as in (15) thru (20) with 

■q replacing \J/ 

/3 replacing 6 



in the matrices, so that 



and 



x m = (CD)„X 



X = (DC),x, 



(26) 



(27) 



where the subscript rj on the matrices indicates the above substitutions. 
Fig. 13 shows the orientation of the earth-centered inclined dipole 
which produces a field of best fit 7 ' 8 to a field based on all observations 
of field vectors made anywhere over the earth's surface. The anomalies, 
or differences between the actual and dipole field, decrease more rapidly 
with increasing heights above the earth's surface than does the dipole 
field itself, making the fit better and better as altitude increases. 
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Fig. 13 — The earth-centered inclined dipole. 



4.5 The Mean Gravity Torque 

If the potential energy of an orbiting satellite is expanded about its 
center of mass in a Taylor series and differentiated with respect to the 
angles giving its orientation with respect to the IS coordinate system, 
the instantaneous gravitational torque acting on the satellite may be 
expressed as 



•here 



T„ = 3 M 2 (/3-/)(vk)(i 8 X k) 



(m 2 = GM/r* 



(28) 



(29) 



GM = universal gravitational constant X the mass of the earth 
r = geocentric distance of satellite 
1 3 = moment of inertia about satellite's spin axis, z 
I = moment of inertia transverse to z (assuming all transverse 

moments to be equal, i.e., I x = I y ) 
k = unit vector along satellite spin axis, z (more generally the 
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vector along the axis of cylindrical mass symmetry, which is 
the spin axis by assumptions in Section 4.1). 

To determine the mean gravity torque, T„ will be integrated over a 
nodal period (7'„) assuming no precession during that period; i.e., the 
SANOR system stays fixed with respect to IS during the integration. 
Then 

(T„) m , UI1 - i- [ ".V(/; t - /i)(vk)(i s X k) dt. (.SO) 

By interchanging the order of the dot product and taking k constant, 

(To) _ = ^^)[ k .{/;!jL. rf( }xk]. (3D 

If a pure Keplerian orbit 6 is assumed, we have 

o(l - c 2 ) 



r = 



r 2 ci = 



GM = 



1 + e cos (<*> — P) 
2ra\l - e 2 ) 1 ' 



T„ 



, 2 3 

V 2 



(32) 
(33) 

(34) 



where 



r = radius vector from the focus of the elliptical orbit to the satellite 
a = semimajor axis of ellipse 

e = eccentricity of the ellipse 
P = argument of perigee 
oj = argument of the satellite. 



Hence the integral becomes 



Jo /•' , Jo 



-— aw 



(35) 



and 
(T.) 



for (A/) 



r„ 2 (i - c 2 ) 3 

x k- / i»i»0 + e cos a> cos P -\- e sin to sin P) dw X k 
Since in the ORDEF coordinates (see Fig. 11) 



(36) 
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i s = cos coi + sin oij g (37) 

the dyad, u* becomes 

i 8 i g = cos 2 oj LA + sin 2 o> jj„ + sin w cos «(ij ff + hh)- (38) 

Substituting this into (3G) results in three integrals. Those containing 
i„i„ and j„j„ yield w as a result of the integration. The other integral is 
0. Therefore, 

(T ) incan = ^^ *' CM. + U) X k]. (39) 

By some vector maneuvering (see Appendix A) this reduces to 

(r.)— . - - «.^ (A7) ,,. [(k-k c )(k X k)]. (40) 

To transform the ORDEF vectors into the SANOR system, we write 

x„ = {CD)& (22) 

X = (DC)x (20) 

.-. x = (CD)(DC)x. (41) 

Performing the indicated operations, after a bit of labor we find that 

kg = (sin ft sin i cos \p — cos ft sin i sin yp)\ 

+ ( — sin ft sin i sin ^ cos — cos ft sin i cos \J/ cos + cos i sin 0) j (42) 

+ (sin ft sin t sin \f/ sin + cos ft sin i cos i// sin 6 + cos i cos 0)k 

where, referring to Figs. 10 and 11, we see that 

yp = an Euler rotational angle 

i = orbital inclination 

ft = ascending node of orbit. 

Substituting into (40), the mean gravity torque reduces to 

(\ 2 /\T 
(T„) mea „ = T2 J _ e8 w ( cos * cos e + sin * sin e cos (0-*)j 

• | (sin % cos cos (ft — $) — cos i sin 0)i + (sin i sin (ft - \p) )]}. 



4.6 The Mean Magnetic Torque 

It is well known 11 that the scalar potential ($>,„) of a magnetic dipole 
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may be expressed as 

•b = — ^ 

4tmo r 2 

-+7TMo r 

= -A^ (45) 

where 

Mu = permeability of free space 

= 47r X 10~~ webers/ampere-meter (Ref. 12) 

= 4tt X 10" 7 henry/meter (Ref. 12) 
p = magnetic moment of earth's field, direction and magnitude 

= 10 weber-meters 

= 10 I7 (10 10 /4ir) emu = 8.0G X 10 25 emu 

= 8.0(i X 10" erg/gauss (Ref. 12, p. 25) 
i„ = unit vector along the x s axis, local satellite vertical 
r = distance from dipole center (earth's center) to satellite 
k,„ = unit vector describing direction of geomagnetic moment 
p = magnitude of p 
r = vector form of r along i a axis. 

A satellite magnetic moment (M) in the earth's magnetic field (H) 
will produce a torque 

T = M X H. (4G) 

If the satellite spins rapidly, any magnetic moment components perpen- 
dicular to the spin axis will tend to produce torques which average to 0, 
while the component along the spin axis will produce a net torque of 

T„, = M s X H (47) 

= MX X H ergs (48) 

where 

M a = satellite magnetic moment along its spin axis (weber-meters 
or ergs/gauss) 
k = unit vector along satellite spin axis (SANOR system) 
H = geomagnetic field (ampere-turns/meters). 

It is desired to integrate the instantaneous torque, T m , over one 
anomalistic period in order to calculate the mean torque. To do this. 
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we shall express H in the MAG system and eventually convert this 
expression into SAN OR terms to be comparable with M s . 
We begin by setting 

H = -V* m (49) 

from (45) and (49), 

R = _h v r-k m -(.r m i m + y m j OT + ^k m )1 (5Q) 

Mo L ** J 

where 

x m , y,„ , z m = components of r in MAG system 

and 

h = 7>/4tt. (51) 

Therefore, 

H=^V( 2m /r 3 ) (52) 

Mo 

and since 

r - (.rj + )j„: + z,,-)'' 

we have 

H = ^ [Sx m zJ m + Zy m z m \ m + ClzJ - r 2 )kj. (53) 

Rewriting to spotlight r components, using the following relationships 
normalized to r, 

x m /r = r x , y m /r = r u , z„,/r = r, (54) 



one obtains 



H = - -A: [3r»rA» + 3r„r j m + (3r, 2 - l)k m ] (55) 



or 



n = -JL-[k m -z(iMQ (56) 

4irMor 

since 

r = ri.. (57) 
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By the following equations [(24) through (65)], the variables i s and 
k„, will be expressed in forms useful in the ensuing mean magnetic torque 
calculations. 

Let us first express i, in the MAG system. As shown previously, 



and 



therefore 



or 



X = (D(~)„X m 
x t = Ai(DC) v X, 



(24) 
(27) 
(58) 



i 3 = | (cos ft cos to — cos i sin to sin ft) 

• (cos to sin ft + cos i cos ft sin to) (sin to sin i) \ 

cos r/ —sin 7? cos /3 sin tj sin /3 
sin 77 cos 77 cos —cos 77 sin /3 
sin cos 

Expanding, 



X 



(59 ) 



i s = [cos 77! cos ft cos to — cos i sin 03 sin ft) + sin t?(cos to sin ft 
+ cos i cos ft sin to)]i„, + [ — sin 77 cos /3(cos ft cos to 

— cos i sin to sin ft) + cos 77 cos /3(cos to sin ft 

+ cos i cos ft sin to ) + sin /3 sin to sin ?']j„, + [sin 77 sin /3 

(cos ft cos co — cos i sin to sin ft ) — cos 77 sin 8( cos to sin 12 

+ cos i cos ft sin to) + cos sin to sin i]k m . 

This will be used later. 

Also, i„ in the IS system may be written [see (24)] as 



(60) 



i. = A ,1 

i, = [cos ft cos to — cos i sin to sin ft[I 

+ [cos to sin ft + cos i cos ft sin to]J 
+ [sin to sin i\K. 



(61) 



(62) 
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In a similar manner, k m in the IS system is expanded as follows: 

x m = (CD) n X (63) 

or 

k,„ = sin 7] sin I - cos 17 sin J -f cos /3 K (64) 

= sin j3(sin r\ I — cos 17 J) + cos /3 K. (65) 

The problem is now to integrate the instantaneous torque, T ffl , over 
one anomalistic period, T a , to obtain the mean magnetic torque, 

(T m ) me a„=^ T m dt=- T m -dV 



Tj H m T m J H m dv 



J a •>!>,• 



(66) 

rln 

h 



where 



T a = anomalistic period of the satellite 
T m = instantaneous magnetic torque 

t = time 

v = true anomaly of the satellite 

/• = geocentric satellite distance 

h = defined by (68) below. 

Substituting into (66) from (48) and (56) 

(T ) = J^-kX r 2 " *= - 3(i -- kJi ' (i) dv. (67) 

U mi mean 4^/^^ ^ ^ J 

Since 

h m W(l - e 2 )* (68) i3 



and 



where 



o(l - e) 
r = 



1 + e cos v 



(69) 



a = semimajor axis of orbit 
e = orbital eccentricity 
T a = anomalistic period of the satellite 
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(07) becomes 

(T ) = M ' V 

V ■*■ mj mean 



87r 2 Moa 3 (l - e 2 )* 

r+ ,, (70) 

•k X / [k m - 3(i,-k m )i.](l + e cos v)dv 

= 37lk X J [k m - 3(i.-kJU(l + e cos v)dv. (71) 

First we shall evaluate the integral of (71) beginning by substituting 
the km value given in (65), to obtain 

(T, "i' nea " = (sin fi)k X r ' [sin v l - cos ij - 3{v (sin ,1 

— cos t?J) }i,](l + e cos y)d« (72) 

+ (cos j8)k X / ' |K - 3(i,-K)i,}(l + e cos w)<fo 

= 1 + 11. (73) 

Integral II, the simpler of the two, is evaluated in Appendix B. The 
result is given below : 

II = 7T cos j8[{ -2 sin $ + 3 sin 2 i sin 

+ 3 sin i cos i cos cos (£2 — ^)ji (88) 

-f- 3 sin i cos i sin (£2 — ^)j]. 

Now we must evaluate the first integral, I, of (73) but first let us 
state that by examining (72) it is perfectly obvious that the I integral 
does not exist for a noninclined dipole. For this case, it follows from 
(73) with cos = 1 that 

\ Am/ mean = ^ 5 ~rrz 9T1 *-*-' (89) 

8ir 2 n a 3 (l — e 2 )* 

By suitable variable substitution the first integral I of (73) may be 
compressed to yield 

I = (sin 0)k X (I. I - JB - 31 C + 3J-E) (90) 

where 

A = / sin 7?(1 + e cos v)dv (91) 
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B = cos tj(1 + e cos v)dv 



(92) 



C = / U. sin f(l + e cosv)dv (93) 

E = / i,i. <'ost?(1 + e cos v)dv. (94) 

Now i? must be expressed as a function of v in order to evaluate A, 
B, C, and E. For simplification of the integrals we shall also let the 
initial time for the integration be the time when the satellite passes 
through perigee. Then 

7J = 7^0 + Ud (9 ,r V 

= tjo + u> E T a M/2T (96) 

= Tjo + bM (97) 

where 

,70 = initial position of the ascending node of the geomagnetic 
equator at perigee passage with respect to the IS system (see 
Fig. 12) 
u B = angular velocity of the earth in the IS system 

I = time measured from passage of satellite through perigee 
M = mean anomaly of the satellite (radians) 
T a = anomalistic period of the satellite (time units) 
I) = (oj B /2x)T rt = number of turns of earth in time T a . 

Therefore 

sin 77 = sin tj cos (bM) + cos ij sin (bM) (98) 

cos tj = cos r/o cos (bM) - sin tj sin (bM). (99) 

Unfortunately, M is related to v through Kepler's equation as 

M = E - eamE (100) 14 

where 

E = the eccentric anomaly 
e = eccentricity of the orbit 
and 



E = 2 tan 



(rr$ ten i] = 2tMrI ( 4 *""!)• (I01) " 
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Equations (100) and (101) certainly define M as a function of v, but 
in a most complicated manner. It can be shown, however, that a plot of 
v and .1/ versus time normalized to T„ will look like Fig. 14. So as a 
reasonable approximation to M we might consider 



M = v - (x/6) sin v 



(102) 



vherc 



\/b = maximum amplitude of the true anomaly "sine" wave of 
Fig. 14 = 2e. 16 

Note that (102) resembles Kepler's equation with v replacing E. With 
less sophistication we might even let 



M = v. 



( 103 



The only justification here is that we shall be dealing in mean torques 
averaged over an orbital period, and the /• function makes one oscilla- 
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280 



240 



200 




160 



Fig. 11 — Mean and true anomaly comparison. 
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tion about the M function in that time. If (103) is assumed, (98) and 
(99) remain unaltered except that v replaces M. If (102) is used, (98) 
and (99) become 

sin r] = sin tj [cos bv cos (X sin v) + sin bv sin (X sin v)] + 

cos ijojsin bv cos (X sin t>) — cos bv sin (X sin v)] 

= sinr/olA'] + cos^B'] (105). 

where A' and B' are denned by comparing (105) to (104). In like 
fashion, 

cos t\ = cos i? [A'] - sin tjo[B']. (100) 

We may express the sin (X sin v) and cos (X sin v) portions of .4' 
and B' as 

» 8 ■ 8 

v • X sin V /m-\ 

sin (X sin v) = X sin v - — ^ — Uwj 

and 

2 • 2 

cos (X sin v) — 1 = • (108) 

These approximations are reasonably valid for near earth satellites 
not having high eccentricities and avoid Bessel function complications. 
The neglect of higher-order terms in (107) and (108) for e = 0.25 and 
T a - 150 minutes results in errors of less than one part in 10 . 

As case I we shall evaluate the integrals A, B, C, and E using (98) 
and (99) with v replacing M. For case II we shall return to evaluate 
these integrals again using (105), (106), (107), and (108). Case I 
details are given in Appendix C. Case II is outlined by Appendix D. 

Using the expanded II integral of (88) and the evaluated A, B, C, 
and E integrals from Appendices C and D, (73) may now be written as 

(T m ) inett „ = 3n|cos/3k X f Q {k - 3(i.-K)i.}(l + e cos v)dv 

+ sin |3kX [Al - BJ - 31- (i„i„Ci + jJ C 2 



+ (U, + 3,0#»)|J 



where 
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A, B, (\ , C 2 , C 3 , Ei , E 2 , E-a are all defined in Appendix C. 

The following, expressed in the SANOR system, will now be substituted 
into (130) (refer to Fig. 11) 

k X I = cos \}/j + sin ^ cos 0i (131) 

k X ( -J) = cos \f/ cos ei - sin i/o (132) 

Ji u = Binfl (133) 

J-j„ = cos Q cos i (134) 

I-i„ = cos SI (135) 

Ij„ = -sin ft cos* (136) 

k X i„ = -cos e sin (Q - \f/)i + cos (12 - ^)j (137) 

k X jt, = ( —sin i sin e — cos i cos e cos (12 — ^))i 

' • • , , < 138 ) 

— cos i sin (12 — ^)j 

k X K = -sin ei. (139) 

These yield the mean magnetic torque, which is 

(Tat) ...em. = 5-5 77T^ ?Ti 1 ^ cos 0[i { — 2 sin + 3 sin 2 i sin e 

+ 3 sin i cos i cos cos (Q — \j/)\ 
+ j{3 sin i cos i sin (12 — r^j] [ 

+ fc g(i%. l^» — + -*■ 

+ Z?[cos ^ cos 0i — sin i/j] 

+ 3[{£ 1 sin 12 + E 3 eos 12 cos i - d cos 12 

+ C' 3 sin 12 cos ij J —cos e sin ($2 — ^)i 

+ cos (12 — i/Oji 4- J/'A> cos 12 cos i + #3 sin £2 

+ C 2 sin 12 cos i — C a cos 0} { — (sin i sin e 

+ cos i cos e cos (12 — ^))i — cos i sin (fi — i/-) j}] [ 

or, collecting on 9ft cos and 311 sin 0, 

(T w )me,.n = tt91Z cos /3 I i{ — 2 sin + 3 sin 2 i sin 0+3 sin i cosi ( 141) 
• cos cos (12 — \p)\ + j| 3 sin i cosi sin (fi — ^)M 



(140) 
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+ 911 sin /3|i j A sin $ cos + B cos ^ cos 

- 3[cos 6 sin (12 — \f/) (2?i sin ft + 2£ s cos ft cos i 

— (\ cos ft + C 3 sin ft cos i) + {sin i sin 

+ cos i cos cos ( ft — ^) ) ( #2 cos ft cos i + #3 sin ft 
+ C 2 sin ft cos i - C 3 cos ft)] f + jJA cos <A 

- B sin V + 3[cos (ft - </0(#i sin ft 

+ E 3 cos ft cos i — C\ cos ft + C% sin ft cos i) 

— cos i sin (ft — ^)(-Sa cos ft cos i 

+ A' 3 sin ft + C 2 sin ft cos i — C 3 cos ft)] j J. 

Using the magnetic moment program, it has been shown in the case 
of the Telstar satellites that letting M = v produces preccssional results 
that follow the case II approximation over 1,000 orbits to within 0.01°. 
To document this result, the residual magnetic moment used was —0.9 
microweber-meter, spin rate equaled 20 to 10 radians per second, spin 
axis and transverse moments of inertia were 4 slug-feet 2 , orbit perigee 
was set to 4,;">00 miles, and eccentricity ranged from to 0.95. Slight 
changes in the x-ij torques of the order of thousandths of a microfoot 
pound were observed as the principal differences in the case I and II 
approximations for these eccentricity ranges. While these differences 
are negligible for the Telstar I and II satellites, other satellites in suffi- 
ciently lower orbits or having greater residual magnetic moments could 
require the case II approximation (no complete study has been made 
to date to bound the required ranges of the above mentioned variables 
for case I to achieve agreement to within 0.01° of case II). 

4.7 Equations of Motion of a Body Symmetrical about Its Spin Axis 

The mean gravity and magnetic torque equations have been derived 
in the previous section. To analyze the motion of a spin-stabilized 
satellite responding to these torques, certain gyroscopic equations must 
now be developed. We begin by relating the vector angular momentum, 
3C, for any rotating body to the external forces acting on the body as 

aC = N (142) 

where 

N = resultant moment of all external forces acting on the body 
3C' = the time derivative of 3C. 
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Using the assumptions in 4.3, we write in the SAXOR system 

3C = ?hk (14.3) 

where 

<p = the satellite spin rate. 
From this, it follows that (see Fig. 10) 

3C' = 7 8 |>k + ?(<» X k)] (144) 

where 

oj = angular velocity of the SANOR system referenced to IS co- 
ordinates. 
We may express o> in the SAXOR system as (sec Fig. 10) 

<d = avi + w„j + w ; k (145) 

or 

<o = 0i + </• sin 0j + ^ cos 0k (140) 

where 

Wj- , u y , a), = angular velocity about the a, //, and z axes respectively. 
Combining (142), (144), and (140), after performing the indicated 
operations 

N = I 3 (frl, sin 0i - <pdj + £k). (140) 

Expressing (140) as x, y, and z torques in the SAXOR system, 

T x = lysine (147) 

T u = -I,tf (148) 

T : = I#. (149) 

Quite obviously it is the x and y torques which produce precession. 
These torques, by the assumptions of Section IV, are the sums of the 
x and I) components of the gravity and magnetic torques expressed by 
(43) and ( 141). 7\ is zero as a result of assuming zero magnetic moment 
transverse to the spin axis and mass symmetry about that axis. We 
note that a nonzero T- implies a change in the spin rate. For the Telstar 
I and II satellites this takes place principally because of induced eddy 
currents which produce transverse moments.* 



* The general equations of motion referred to the center of mass for a rigid 
body .spinning about the z axis, and symmetrical about this axis, are 
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4.8 Alternate Inertial Coordinate System 

Clearly, (143) is valid only if the angular momentum vector, DC, 
coincides with the satellite spin axis k. We note in (147), which 
was derived from (143), that the above assumption will be approxi- 
mated only if \j/ sin is small as compared to <p. Otherwise the direction 
of 3C shall certainly be influenced by that term as well as <p. We note 
too that a valid condition within the bounds of the assumption is for 
^ sin to be small even with large ^ provided only that itself be ap- 
propriately small. That is, solutions to the equation of motion exist 
for large \f/, and those will occur only for small because of (143).* 
But if be near 0° or 180°, a singularity exists in (147), for \j/, as a 
result, approaches infinity. An exit from this dilemma may be secured 
simply by transforming to a new reference set of inertial coordinates in 
place of the IS system whenever becomes small. Naturally, the new 
set should be chosen so that the equivalent then existing will be 
large. This is accomplished by redefining IS so that 

Xj corresponds to Y 

Y-i corresponds to Z 

Zt corresponds to X 

as shown in Fig. 15. Let us call this new inertial frame the IS2 system. 
Transformation equations to relate IS to IS2 are quite simple and 
are given below : 

X 2 = QX (150) 

X = QX 2 (151) 



T z = Id + (h - IW sin e cos e + lv(4> sin 

T y = H sin + (21 - U)W cos - IvpO 
T z = h{$ cos - <j>6 sin + v>). 

For cases where <p dominates , 0> sin 0, and \p cos 0, these torques reduce to 
those given by (147), (148) and (149). We note principally that a change in spin 
rate reflects a T x processional torque, from the above equations, but this is 
generally small and is neglected in this paper. For Telstar I, this torque compo- 
nent is estimated at least two orders of magnitude below the magnetic torques 
considered herein. (See also Ref. 17.) 

* This is the case for precession through or near the north or south celestial 
pole, where even small changes in attitude, or 0, can produce large changes in ^. 
The situation is quite analogous to an azimuth-elevation antenna tracking a 
satellite that passes through or near the zenith, where the azimuth rates become 
extremely high even for small changes in satellite position. 
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Fig. 15 — Euler rotations in the alpha-gamma system. 



Q = 



and 



Q = 
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(152) 



(153; 



Relationships between the other coordinate systems and the IS2 system 
may be obtained using equations of Section 4.4 and applying the Q 
matrix as appropriate, except that in the expanded matrices 

a corresponds to \J/ 

7 corresponds to 6 

where 

a = The Euler angle measured from X 2 to the intersection of the xy 
and X 2 F 2 plane, which defines a new x axis called the x 2 axis in 
a SANOR 2 system. (See Fig. 15.) 
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7 = The Euler angle measured from Z-< to k in a plane perpendicular 
to the Xi axis. ( Rotation is about the .r 2 axis. ) 

Quite obviously, the steps leading to the mean gravity and mean 
magnetic torques may be retraced using the Q matrix and the new 
a,y Euler angles. This would lead to torque expressions as functions 
of these desired angles. The same result may be obtained by transforming 
the final torques expressed in the yp-6 system to expressions in the a-y 
system by using explicit relationships between these two systems. We 
shall choose this later route and the needed relationships will now be 

derived. 

Wc note in Fig. 15 that the k (satellite spin) vector remains in the 
same spatial position whether expressed in the IS or IS2 system. It 
seems reasonable then to proceed to relate a-y to ifr-O by observing the 
projections of the k vector on the XY and X 2 F 2 planes, respectively. 

Referring to Fig. 10, we sec that 



tan a = 



— sin 6 cos \p 



— cos 



cos 7 = sin sin \f/ 

sin 7 = (1 — sin 2 sin \p) . 



(157) 

(158) 

(159) 
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Fig. 16 — Theta-psi and alpha-gamma geometric relationships. 
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Quite obviously, inverse relationships may he developed. These are 

tan* = — . cos 7 . (1G0) 

— sin a sin 7 

cos 6 = —cos a sin 7 (161) 

sin = (1 - cos 2 a sin 2 7)*. (162) 

•t.9 Mean Gravity and Magnetic Torques — Alternate Euler Rotations 

Using (157) through (159) it can be shown that the mean gravity 
torque expressed in the a-y system is 

/™ \ 67r 2 A/ I , . _ . . 

(T«) mean = jrjn _ 2 )i J C— sin iz sin 1 cos 7 

+ cos ft sin i sin a sin 7 + cos i cos a sin 7) [ 

. . . ... * (163) 

[i 2 (sin ft sin i sin 7 + cos ft sin 1 sin a cos 7 

+ cos i eos a eos 7) + jy(cos ft sin i cos a 

— cos i sin a)] . 
Likewise, the mean magnetic torque can be expressed as 
(T m ) mcB n = 7ri)llcos/3 vA — 2 cos a cos 7 + 3 sin i( sin i cos a cos 7 

— sin ft cos i sin 7 — cos ft cos i sin a cos 7) r 
+ j>j2 sin a — 3 sin i(sin i sin a 

+ cos ft cos z cos a) M 

+ 9TC sin /i I.] —.1 sin 7 

L . (1C4) 

— B sin a COS7 -+■ 3Gi(sin ft sin a cos 7 

— cos ft sin 7 ) + '.\(1<{ —sin £ cos a cos 7 

+ sin ft cos i sin 7 + cos ft cos i sin a sin 7) r 

+ h\—B eos a + 3(?i sin ft cos a 

+ 36' 2 (sin i sin a + eos i cos ft cos a) M 

where 

(V, = Ei sin ft + #3 cos ft cos i - (\ cos ft + C 3 sin ft cos i (165) 

C/j = 7?o cos ft cos i + /'>' :i sin ft + C2 sin ft cos i — (\ cos ft. (166) 

All other terms have meanings described previously. 
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V. EXPERIMENTAL RESULTS — GENERAL 

Using the methods described herein, the attitudes of the Telstar I 
and II satellites have been predicted, determined, and these two results 
compared. In the case of the Telstar I satellite, the maximum angular 
deviation between predicted and determined attitude is 0.90° with the 
special-case exception of pass 2154, which is covered in detail later.* 
The average angular deviation is 0.38°. The predicted attitude of the 
Telstar II satellite deviates from that determined by optical flashes and 
solar sensors by a maximum of 0.5° and an average of 0.09°. These 
results are tabulated in Tables II and III. Here the deviation in right 
ascension and declination is displayed by subtracting the predicted 
values from the determined values. The angular deviation between pre- 
dicted and determined attitude is also given as the celestial great circle 
arc G. When this value is followed by N, this indicates that no attitude 
fix has been determined by loci intersections (see Section II), and so 
the shortest arc distance between the corresponding attitude locus and 
the prediction is quoted. The T preceding an attitude determination indi- 
cates the selection of that point as a target toward which a magnetic 
moment program embodying the above analysis attempts to converge. 

5.1 Experimental Results — The Telstar I Satellite 

The Telstar I satellite entered orbit on July 10, 1962. Its attitude, 
calculated by combining nominal third-stage burnout parameters and 
stage 1 and 2 telemetry data, is given in Table I. Between launch and 
pass 16, there were a number of instances when the orientation coil was 
inadvertently energized by misinterpreted command signals sent to the 
satellite at extreme ranges and/or low elevations. This orientation coil 
(often referred to as a "torque coil") is located in the equatorial plane 
of the Telstar satellite just under the outer skin. Its purpose is to enable 
attitude adjustments by the production of a magnetic moment whenever 
current is sent through the coil. A fixed amount of current may be caused 
to flow in either direction through the coil so as to produce a magnetic 
moment, in addition to the residual moment, of ±7.8540 microweber- 
metcrs. Fig. 17 diagrams the sense of the magnetic vectors which are 
called positive, indicates the resulting north and south magnetic poles 
of the satellite, and shows the corresponding positive direction of the 
current flow in the orientation coil. The relative magnetic directions and 
current flow are in accord with established standards. 



* There are apparent exceptions on passes 16, 72 and 1657, but these relate to 
antenna pattern techniques (see Section 5.1) which have an expected accuracy 
of ±1° to ± 2°. 
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Table II 


— Telstar I Satellite — Attitude Comparison 






Attitude Deter- 
mination 


Type 

Data 


Attitude Prediction 


Deviation 
(Det. — Predict.) 


MAG 


Pass Number 




for 
De- 






MOM 
















(MWM) 




RA 


DEC 


term. 


RA 


DEC 


RA 


DEC 


G 




Burnout 






83.7 


-00.8 










10 


84.2 


-05.8 


SM 


84.2 


-05.18 











-0.80328 


16 


83.0 


-00. 


A 


84.2 


-05.18 


-1.2 


-0.82 


1.4 






72 








87.78 


-60.27 






0.5N 






72 


80.0 


-67 


A 






-1.3 


-0.73 


1.2 






135-130 


91.7 


-00.0 


M 


92.19 


-06. 3G 


-.49 


0.36 


0.0 






199 


single locus 


M 


90.49 


-05.95 






0.4N 






271 


100 


-05 


A 


100.80 


-04.92 


-0.8 


-0.8 


0.5 






272 


single 


OCUfl 


M 


100.90 


-04.90 






0.0N 






272 


I 


I 


S 


100.90 


-04.90 






0.2N 






472 




107.50 


-59.80 






0.2N 






931 


T100.5 


-49.84 


M 


100.74 


-49.93 


-0.24 


0.09 


0.18 
























-0.80328 


1051 


single locus 






97.25 


-50.21 






0.9N 


+0.9907 


























-0.80328 


1009 














98.42 


-49.70 






0.GN 


-8.7173 


1114 














95.92 


-51.01 






0.1N 


-0.71103 


1430 














92.02 


-59.55 






0.1N 


I 


1507 


T96.2 


-64.2 


M 


90.71 


-04.21 


-0.51 


-0.01 


0.23 




















-0.71103 


1G57 


99.2 


-67.0 


A 


102.81 


-66. 64 


-3.01 


0.30 


1.7 


7.14297 
-0.59101 


1095 








78.58 


-57.01 












1909 








90.59 


-59.71 






0.7N 






2154-55 


T104.0 


-53.5 


M 


104.71 


-54.81 


-0.71 


-1.3 


1.5 






2200 


single locus 






100.03 


-53.28 






0.1N 






2204 














107.20 


-50.74 






0.4N 






2404 














107.09 


-43.05 






N 






2482 














100.88 


-43.04 






N 






2509 














100.57 


-42.33 






N 






2582-83 


105.2 


-40.2 


M 


105.44 


-40.53 


-0.24 


0.38 


0.4 






3340, 41, 42 


100.3 


-43.7 


M 


103.08 


-01.54 


-2.02 


17.8 


18.0 






3470-7 


112.30 


-40.32 


M 


120.21 


-00.21 


7.91 


19.9 


22.0 






3495 








123.10 


-00.41 












3470-7 


T112.30 


-40.32 


M 


112.30 


-40.32f 











-0.40945 



M = Mirror (lash data. 

A = Antenna pattern data. 

8 = Solar sensor data, 
t Connects 3340 to 3470. 



Data relating to the inadvertent uses of the orientation coil are insuffi- 
cient to reconstruct the detailed processional motion of Telstar I from 
launch to pass 16. Cn pass Hi, however, an attitude determination was 
made by reducing optical and solar sensor data. This along with an 
attitude fix on pass 135-136 made from mirror flash data indicated a 
magnetic moment of about —0.7 microwcber-meter. This was later 
refined to — 0.8G328 microweber-meter by causing the magnetic moment 
program to connect, by the processional theory herein described, the 
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Table III— Telstar II Satellite — 


Attitude Comparison 




Puss Number 


Attitude Determin- 
ation 


Type 
Data 


Attitude 
Prediction 


Deviation 
(Dct.-Predict.) 


MAG 
MOM 

(MWM) 




RA 


DEC 


RA 


DEC 


ARA 


ADEL 


G 


Burnout 
62, 63, 63 


87.75 


-55.08 


M 


82.23 

87.75 


-57.31 
-55.08 











-0.48375 


100 

132, 133, 133 


88.75 


-55.34 


M 


88.55 
89.23 


-55.05 
-54.97 



-0.48 



-0.37 



-0.5 


-0.48375 


281 

293, 293 
324, 25 


92.3 
92.40 


-54.05 
-53.71 


M 

M 


92.17 
92.38 
92.90 


-54.10 

-53.99 
-53.71 


-0.08 



-0.0G 



ON 
0.1 








331 








93.01 


-53.66 






0.08N 






490 
496 
541 
573 
611 
643 


T95.00 


-51.73 


M 


94.95 
95.00 
95.34 
95.49 
95.63 
95.71 


-51.81 
-51.73 
-51.10 
-50.74 
-50.20 
-49.80 






ON 

ON 

ON 

0.13N 
0.16N 


> ' 



M = Mirror flash data. 
S = Solar sensor data. 

attitude on pass 16 to a point on the pass 931 attitude locus. The con- 
vergence is within 0.18°, as shown in Tahle II. Using this latter magnetic 
moment, the precessional motion from passes 10 to 1051 was established. 
This precessional history is given in Fig. 18, where both the predicted 
curve and the attitude loci data are displayed. 

When an optical flash series is recorded at the Holmdel Laboratories, 
the midpoint of that series can be determined with a time accuracy of less 
than ±20 seconds when using visual observation through the on-site 
telescopes and within ±10 seconds photoelectrically. The early data 
recorded on the Telstar I satellite did not indicate the expected accuracy 
at the midpoints, so in Fig. 18 only central loci are plotted. We must 
therefore think of these loci as having possible tolerances up to ±10 
seconds, since photoelectric data were reduced. A tolerance of ±10 
seconds can cause the maximum locus limit to be displaced as much as 
±0.7° on either side of the plotted central loci. Considering this, the pre- 
dicted attitude curve passes through every loci and every loci intersec- 
tion on Fig. 18 with the single exception of pass 1051. The prediction 
curve would miss the lower limit of this locus by 0.3°. 

Extended attitude predictions were made on October 1, 1902 (around 
pass 800) and covered the period through January 22 (pass 1800). 
These indicated that attitude adjustments should be initiated during 
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January, 1963. This was needed to prevent the solar offset from exceed- 
ing about 15° as required (see Section I). In order to check the attitude 
correction procedures, a test maneuver was initiated on pass 1051. It 
began by turning on the orientation coil, so as to produce a magnetic 
moment of +7.8540 microweber-meters to oppose the -0.86828 residual 
moment and yield a resultant of +6.9907 microweber-meters. This was 
continued through pass 1058, whereupon the coil was turned off until 
pass 1069. On 1069 the coil was turned on in the opposite sense so as to 
produce a moment of —7.8540 microweber-meters, giving the satellite 
a net moment of -8.7173 microweber-meters. On pass 1075, the coil 
was turned off. All told, the coil was positive for about 18 hours and 
negative for about 16 hours, a fact which approximately nullified the 
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Fig. 17 — Direction of positive magnetic moments for Telstar satellites. 
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Fig. 18 — Telstar I precessional history — passes 16 through 1051. 

precessional effect of the coil as intended for this test. Details of this 
maneuver are shown in Fig. 19. 

Table II summarizes the agreement of loci to prediction during this 
period. Giving loci of unknown tolerances a ±0.7° spread and using the 
known tolerance for pass 1114, the residuals shown in Table IV are 
produced. 

The magnetic moment program connected the attitude on pass 1075 
to that determined by mirror data on pass 1567 within 0.23°. In so 
doing, it calculated a residual magnetic moment of -0.71103 micro- 
weber-meter for that era. Using this calculated moment, predictions are 
extended to pass 1657 where, on January 7, 1963, the main orientation 
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TO PASS 1567 
WITHIN 0.23° 
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Fig. 19 — Details of the trial torqueing maneuver. 



Table IV — Telstar I Satellite — Test Maneuver 



Pass Number 


Deviation of Prediction from Optical Locus 


1051 
1069 
1114 


0° 
0.1° 
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maneuver began. The region from pass 1075 to 1657 is shown in Fig. 20 
along with the attitude loci. Predictions fit the determinations through- 
out within about 0.2°. 

In addition to using mirror flash and solar sensor data to determine 
attitude, W. C. Jakes, Jr. and group at the Holmdel Laboratories 
deduced the attitude for a number of passes in this period by analyzing 
antenna pattern data. 18 Some of these determinations are included in 
Table II. Since their expected accuracy is ±1°, these data also serve 
to corroborate the predicted attitude. 

The orientation coil was energized for the major attitude maneuver 
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Fig. 20 — Telstar I precessional history — passes 1075 to 1657. 
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beginning on pass 1657, January 7, 1963, and ending on pass 1695, 
January 11, 1963. During this time the residual magnetic moment of 
— 0.71103 microweber-meter had superimposed upon it a 7.8540 micro- 
weber-meter moment as a result of the coil. The attitude of the Tclstar 
I satellite was changed by over 20° during this maneuver. The attitude 
and curve for this period are given in Fig. 21. The slight ripples in these 
curves are due to the inclination of the earth's dipole. Because of the 
small scale of this plot, these ripples are more noticeable than in the 
previous graphs. 

At this point, we might well insert a brief discussion of the planning 
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Fig. 21 — Telstar I precessional history — passes 1657 to 1695. 
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behind this attitude maneuver and why pass 1657 was chosen for its 
commencement. First of all, it can be shown that for any angular offset 
in the spin axis from a fixed, desired orientation there exists, for each 
instant in time, a determinable polarity of voltage for the equatorial 
coil which will tend to decrease the offset. To utilize this principle fully 
would require turning on the orientation coil on many occasions when 
the Telstar I satellite entered the Andover skies. While such a procedure 
would tend to decrease the angular offset generally in the most direct 
fashion, there are occasions in the lifetime of the Telstar satellite when 
leaving the coil on for a number of revolutions will produce faster cor- 
rections than pulsing the coil continually at Andover. This is because the 
instantaneous torque produced by the continuous magnetic moment of 
the coil in the earth's magnetic field varies during an orbital period in 
such a way that, in general, the average value exceeds the torque obtain- 
able by energizing the coil only in the Andover skies. The optimum pro- 
cedure then is for a period of continuous operation and this was under- 
taken. Fig. 22 shows the results of turning on the coil for steady torqueing 
at various dates. The heavy spiral shows the predicted position of the 
spin axis in right ascension and declination if the coil is never energized 
(the trial maneuver is omitted here for simplicity). The diverging curves 
indicate the motion due to the torqueing coil. Examination of Fig. 22 
shows rapid precession as a result of the torque coil field, so that no mode 
of continuous coil operation will long permit the spin axis to point close 
to the south ecliptic pole. The solar offset also must ultimately increase 
with the enlarging ecliptic angle. It therefore follows that if a simplified 
mode of coil operation exists, it must begin with a relatively short period 
of coil operation. 

Two factors concerning the precessional motions are noted in Fig. 
22. They are: 

(i) All precessional motions involving the orientation coil lie outside 
the spiral of residual precession for many orbits, because the net mag- 
netic moment existing when the coil is actuated exceeds the residual 
magnetic moment. 

(ii) Torque coil precession for both positive and negative coil polarity 
begins incrementally at the spiral in opposite directions. 

We see that the torque coil should be used sparingly, because of the 
relatively rapid motion it produces. The problem that remains is to 
determine the times to turn on and turn off this coil and also to investi- 
gate whether or not future coil use is required. Because the spin rate of 
the satellite is decreasing, equal magnitudes of satellite magnetic mo- 
ments will cause greater and greater precessional motion as time passes. 
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Fig. 22 — Telstar I predicted attitude showing stead}' torqueing. 



It is to be expected that after the coil is actuated and then turned off 
a spiral of larger excursion than that occurring for passes 16 through 1655 
will result strictly because of the lower spin rate existing at the later 
time. There is no way to avoid this except by repeated coil pulsing. 

To avoid, or at least minimize, repeated coil usage, the procedure is 
to try to keep the maximum ecliptic angle as small as possible. If this 
value is expected to exceed the allowable solar offset, it is advantageous 
to time the operation so that the maximum ecliptic angle will occur at 
the time of minimum solar offset. This criterion governs both the selec- 
tion of the day for coil turn-on as well as the duration of the attitude 
maneuver. As may be guessed from Fig. 22 and from the processional 
motion factors stated above, only short torqueing operations during the 
month of January would result in the next precessional spiral (after coil 
turn-off) returning anywhere near the initial spiral shown in Fig. 22. 
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To cause the maximum ecliptic angle and the minimum solar offset to 
occur together, orientation coil use had to begin in the period from 
January 1 to January 7, 1963. As shown above, the maneuver which 
did begin on January 7 placed the attitude at right ascension 78.6° and 
declination at 57.0° by pass 1695. 

The magnetic moment program, in connecting pass 1695 attitude to 
that determined by mirror data on passes 2154 and 2155, determined a 
residual magnetic moment of —0.59101 microweber-meter. This differs 
by about 0.1 microweber-meter from the residual moment going into the 
torqueing maneuver, but appeal's to be borne out by the closeness of 
fit to the mirror data from pass 1695 through pass 2583 (see Table II).* 
Fig. 23 plots this trajectory along with the mirror loci. The fit of attitude 
predictions to determinations is very good from passes 1695 through 
2583 (see Table II). The largest deviation of 1.5° on passes 2154-2155 
may be somewhat misleading. Referring to Fig. 23, we see that the 
predicted attitude on pass 2154 is easily within 0.1° of the 2154 locus. 
The interaction of the 2154 and 2155 loci however, produces a common 
area, or attitude box, 1.5° from the pass 2154 prediction. This box is 
produced by a somewhat grazing intersection of the loci, and its location 
is therefore affected by the precession which took place between these 
two passes. Antedating the 2155 locus to account for this would place 
the attitude box within 0.05° of the 2154 and 2155 prediction, so this 
fit is quite valid. 

We are not in so comfortable a position for passes 3340 and there- 
after. Attitude boxes for 3340-1 and 3476-7 arc shown as dashed lines 
on Fig. 23. They are about 20° from the corresponding predicted atti- 
tudes, and the reason for this discrepancy is at present unknown. We 
note that the Telstar I satellite ceased its transmission on pass 2065 
(February 20, 1963) because of radiation damage, but that the attitude 
remained predictable at least through pass 2583. The region from pass 
2583 to 3340 is devoid of attitude data because no telemetry could be 
received to report solar aspect and no mirror flashes were recorded. This 
was due to the increased activity on Telstar II, a certain proportion of 
Telstar I passes occurring in daylight hours, and prevailing weather 
conditions at Holmdel. 

Interestingly enough, it is noted that the attitude determinations on 
passes 3340-1 and 3476-7 can be connected by the magnetic moment 
program even though they cannot be sensibly joined to the preceding 



* It is possible that coil usage can alter the residual moment of the satellite, 
which was made small in the first place by appropriately balancing much larger 
magnetic fields. 19 
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60. 90. 

SPIN AXIS RIGHT ASCENSION IN DEGREES 



TQ.STAR I ATTITUCE - BELL TELEPHCfC UBOMTCRIES - L C TKMAS 
MAG MX = -.59101 »**, CRBIT ASI2 



120. 



Fig. 23 — Telstar I prer-essional history — passes 1G95 to 3500. 



data. Fig. 24 shows this union. The magnetic moment of —0.59101 
microweber-meter used in these last three plots cannot be considered too 
reliable, since it is one of a number which can connect the two attitude 
boxes shown in Fig. 24. However, it is the same as previously used in the 
1695 to 2583 region and may tend to be valid because of this. In any 
event, more attitude data are required to resolve this issue. 

As to what may have caused the anomalous behavior of Telstar I 
somewhere between pass 2583 and 3440, we offer the following, taken 
singularly or in combination, as possibilities: 

(i) the orientation coil has been energized during this period (from 
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Fig. 24 — Telstar I precessional history and prediction — passes 3340 to 5000. 



the present precessional motion, it appears to have been turned off 
between 3340 and 3477) ; 

(it) the residual magnetic moment has changed a multiplicity of 
times; 

(Hi) meteoric collision has taken place; 

(iv) pressure leakage from instrument canister has occurred in a 
manner to alter the attitude by reaction forces. 

5.2 Experimental Results — The Telstar II Satellite 

The Telstar II satellite entered orbit on May 7, 1963. Its initial 
attitude calculated from the third-stage burnout parameters is shown 
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in Table III. The predicted burnout attitude proved inaccurate because 
the third stage of the Thor-Delta vehicle did not perform nominally. 
This is evidenced by the fact that the orbital period calculated from 
such assumed nominal performance differed from the actual period by 
as much as four minutes. 

The first attitude fix of the Telstar II satellite occurred on passes 62 
and 63, when a total of three flash series were observed. Connecting 
this determination to that of pass 496 requires a residual magnetic 
moment of —0.48.375 microweber-meter. The fit through pass 643 is 
given in Table III, and some typical attitude loci arc plotted along with 
the attitude prediction in Fig. 25. Corresponding ecliptic angle, solar 
aspect, X and Y torques appear in Figs. 26 through 29. 




e». 90. 92. 94. 96. 90. 

SPIN AXIS RIGHT ASCENSICN IN DEGREES 



TELSTAR II ATTITUDE - BELL IELEPMCTC LABCHATCRIES - L C THCHAS 
MAGNETIC MrMENT=-.A8375 H-W-M, OT13, INITIAL ATT FM FLASHES PASS 62,63 



Fig. 25 — Telstar II processional history and prediction — passes 62 to 3000. (See foot- 
note, p. 1713.) 
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PERIGEE NIWER 

TELSTAR II ATTITUDE - BELL TELEPHCtC LABORATORIES - L C THCMAS 
MAGNETIC *>CNT=-. 48375 K-1/-M, OT13, INITIAL ATT FM FLASHES PASS 62,63 

Fig. 2G — Telstar II ecliptic angle — passes 02 to 3000. 



The solar offset is due to exceed the 15° limit dictated by temperature 
balance considerations (see Section I) by pass 1250 (see Fig. 27). Since 
the maximum offset which occurs on pass 1400 is only 18°, no plans are 
contemplated to reorient the Telstar II satellite strictly to prevent this 
mild excursion. 

There is, however, another excursion beyond the 15° limit, around 
pass 2670. Besides this excursion being more serious than the former, 
some interesting differences between these two events exist, as Fig. 
25 illustrates. 

On pass 1400, the attitude will be returning to a region nearer the 
south ecliptic pole (located at right ascension 90°, declination 67.5°) 
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TELSTAR II ATTITUDE - BELL TELEPtOC LAKSATCTJIES - L C TKMAS 
WWCTIC M>CNT=-. 48375 M-lrf-M, OTI 3 • INITIAL ATT FH FLASHES PASS 62,63 



Fig. 27 — Telstsir II solar aspect — passes 62 to 3000. 



thereby Limiting the solar offset to values below 15°. On pass 2670, 
however, the attitude is moving away from the south ecliptic pole, and, 
in fact, because of its now lower spin rate is entering a spiral of greater 
excursion than that previously traversed. It therefore behooves us to 
correct the attitude before that latter spiral occurs. 

The question remaining is to determine the most profitable time for 
such a correction. Ideally, a well-chosen time would meet the following 
conditions: 

(i) It would take place just after a good optical attitude fix was 
established to verify the predictions. 
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Fig. 28 — Telstar II A' torques — passes G2 to 3000. 

(ii) The attitude maneuver would result in as many mirror flashes 
as possible during correction to indicate its progress. 

(Hi) The attitude maneuver would end with the satellite in a posi- 
tion to guarantee a suitable attitude for as long a time in the future as 
is possible. 

At the time of this writing, suitable attitude corrections are under 
study. Most probably the orientation coil will be energized sometime 
between passes 1300 to 1500 in a negative sense or in a positive sense 
near pass 2400 so as to drive the attitude downward and to the left in 
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Fig. 29 — Telstar II Y torques — passes 02 to 3000. 

the sense shown in Fig. 25. This permits the next attitude spiral to occur 
near the south ecliptic pole.* 



VI. CONTUSIONS 



Techniques for both attitude determination and prediction for spin- 
stabilized satellites have been developed. Their use has been demon- 

* The attitude of the Telstar II satellite was successfully reoriented by energiz- 
ing the coil on pass 2402 (May 17, 1904) and leaving it in that state until pass 2421. 
Thus attitude predictions given in Fig. 19 beyond pass 2402 are no longer valid. 
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strated using Telstar I and II satellite data. It has been shown that an 
inclined dipole model of the earth's magnetic field and the method of 
averaging the gravitational and magnetic torques over each anomalistic 
period of the satellite permit attitude predictions to within a few tenths 
of a degree of determined values in most instances. In those few cases 
where departures are above one degree, explanations have been pre- 
sented to show the reason for such discrepancies. The reasons are (1) 
unknown time errors in determining the midpoint of the optical flash 
series and (2) grazing intersections of the attitude determining loci. 
There remains but one anomaly in the precessional motion of the Telstar 
I satellite which, for the moment, is unexplained. Possible reasons for 
this anomaly are given in Section 5.1. 

It has further been shown that the seemingly crude approximation of 
letting the mean anomaly of the satellite equal the true anomaly for the 
purposes of determining the mean torques produces attitudes in close 
agreement (0.01° over 1000 orbits) with more sophisticated approxima- 
tions for orbit eccentricities up to 0.9 and perigee radii above about 
4500 miles. (See Section 4.6, case II.) 

The comparisons made herein of the precession given by the magnetic 
moment program and the attitude determinations substantiate the 
simplifying assumptions made in Sections 4.1, 4.2, and 4.3. These 
assumptions are most instrumental in producing a working technique 
for both attitude prediction and residual magnetic moment determina- 
tion which is amenable to analytic solution and conservative of com- 
puter time. 

Furthermore, the usefulness of combining optical flash and solar 
sensor data for attitude determination and their inherent accuracy is 
shown. Optical flash data can provide loci with a resolution of 0.1°. 
Solar sensor loci are resolved to within 1°. While it is clearly straight- 
forward to determine analytically the boundaries of the attitude boxes 
from intersecting loci exhibiting estimated time tolerances, this paper 
indicates the decided advantage of graphing the individual loci to 
determine the angle of intersection and thereby gain an estimate of the 
validity of the boundaries in the presence of precession (see, for ex- 
ample, Section 4.1). 

Finally, the techniques described have all been consolidated into 
working computer programs which follow closely the analysis presented. 
In addition, a number of important supporting calculations such as the 
solar position, sidereal time, orbit updating, etc. are developed. Because 
of the complexities of the mean torque and gyroscopic equations, the 
precessional analysis is most useful when embodied in suitable computer 
programs. 
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APPENDIX A 

Derivation of Equation (40)" 

Beginning with (39) of the main body of the paper, it is certainly 
evident that 

k-(i tf i„ + j„j ) X k = (k-i„)(i„ X k) + (k-j (/ Hj„ X k). (39a) 
It will be shown that the following identity exists 

(k-i.Xi, Xk) + (k-j„)(j fl X k) + (k-k„)(k Xk)=0. (39b) 
Rewriting (39b), one obtains 

[(k-iji, + (k-j„)j fl + (k-kjkj X k = (39c) 

and 

(k-i„) = (:r„i + y o j + z \s„)-i B (39d) 

= x„ (39e) 

where 

x a ,y<j j z a = components of k along the ORDEF axes. 
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Iii like fashion 

(k-J.) = 'Jo (39f) 

(k-k a ) = z a . (39g) 

Substituting (39e,f,g) into (39c), we obtain 

(x„i g + y„h + zX) XkEO (39h) 

k X k = (39i) 

= 0. (39j) 

Therefore 

-(k-k„)(k fl X k) m k-(Mo + hh) X k, (39k) 

and (40) follows from (39). 

APPENDIX B 

Evaluation of the II Integral of Equation (73) 2> 

Copying the II integral from (73) of the main body of the paper, we 
have 

II = (cos j8)k X ( (K - 3(i.-K)i.) (1 + e cos v) dv. (74) 
Split (74) into two parts by letting 

IIA = cos 13k X K [ (1 + cos v) dv 

Jo 

= cos/3K X / (1 + e cos a> cos P + e sin a> sin P) dco 

Jo 

= (27rcosj3)k X K 
and 

IIB = (-3cos/3)kXK- / i.i.(l + e cos « cos P 

J o (76) 

+ e sin to sin P) dco . 

But 

i„ = cos wi„ + sin wj„ (77) 
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therefore 

i s i s = cos 2 tAgig + sin 2 «jj„ + sin co cos w(i„j„ + j ff i„). (78) 

It follows that 

IIB = (-3tt cos/3 )k X K-(i„i fl + jJJ. (79) 

So that, combining IIA and IIB 

II = (7rcos/3)k X [2K - 3K- (i„i„ + j„j„)] (80) 

= (xcos^)k X [2K - 3(K-j )j„] (81) 

= 7rcos/3[2k X K - 3(K.j„)(kZj p )]. (82) 

We shall now proceed to express (82) in the SANOR system. From 
Fig. 10, 

K = sin 0j + cos 0k. (83) 

From (22) and (20) we write 

x a = (CD)i(DC)x (84) 

or 

j B = (sin yp cos 12 cos i — sin 12 cos \// cos i)'\ 

+ (sin i sin + sin £2 cos i sin ^ cos 

+ cos 12 cos i cos i/< cos 0)j (85) 

+ (sin i cos — sin 12 cos i sin \f/ sin 

— cos 12 cos i cos ^ sin 0)k 
or, simplifying 

j„ = [— cos i sin (12 — ^)]i 

+ [sin i sin 4- cos i cos cos (12 — ^)]j (86) 

+ [sin i cos — cos i sin cos (12 — ^)]k. 
Substitute (83) and (86) into (82) to obtain 

II = 7T cos /3[( -2 sin 0)i - 3[jsin i sin 2 

+ cos i sin cos cos (12 — xp) + sin ?' cos 2 

— cos i sin cos cos (12 — $)} { ( —sin i sin (87) 

— cos i cos cos (12 — \p))i 

+ (-cost sin (12 - i/0)j)]]. 
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Finally, 

II = 7T cos 0[{ -2 sin + 3 sin 2 i sin 6 

+ 3 sin » cos i cos cos (8 — ^)}i (88) 

+ 3 sin i cos i sin (fl — ^)j]. 

APPENDIX C 

Case I Expansion 

Refer to integrals of (91), (92), (93), and (94) of the main body of 
the paper: 

.2* r 2 "- r 2 * 

A = I sin 7?o(l + e cos v) dv = / sinijodo + e / sin ijo cos » dv 

'' " ° (109) 



= sin 770 

'0 



/ cos (bv) dv + cos 770 / sin (bv) 
.'0 J o 

■f e sin i/o / cos (&v) cos v dv 

1° (U0) 

+ cos 170 / sill (bv) cos v dv 

= sin tfi(Di + eD 3 ) + cos 770 (A> + eD 4 ) (HI) 

= ( I f^_ } [ s i n ^o s in 2ir6 + cos 770 (1 - cos 2tt6)] (112) 

\b 1 — o 2 / 

where the D factors are listed in Appendix E. In similar maimer, 

B = cos ij (Di + eD 3 ) - sin ij (D 2 + eD 4 ) (113) 

= H ^L_l [cos 770 sin 2tt6 - sin 770(1 - cos 2tt6)] . (114) 

The C integral is expanded as follows 

f.2* 

C = / U s sin 77(1 + e cos v) dv (115) 

Jo 

but 

i„ = cos wig + sin &>j„ . (116) 
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Therefore 

y, = cos 2 wi„i B + sin 2 wjja + sin w cos u ^» + ^' i|; ^ (U7 ^ 

Let 

c « ma + ua + cu, + y.)Ci < 118 ) 

and evaluate C, , C 2 , C a separately. 

f\ = /" cos 2 u sin ij(l + e cos v) dv . (119) 

Jo 

Letting co = t> + P, expanding cos 2 co and collecting terms, we have 
C, = sin Wcos 2 P(D 8 + eD 7 ) + sin 2 P(D 9 + eDn) 

- 2 sin P cos P(Di 3 + eD n )] , 120 ^ 
+ cos r;olcos 2 P(D 6 + rD 8 ) + sin 2 P(D* + eDu) 

- 2 sin P cos P(D M + eDu)]. 
in like fashion 

C 2 = sin rjolsin 2 P(D 5 + eD 7 ) + cos 2 P(A + eDu) 

+ 2 sin P cos P(Du + eDi B )] , 121 ) 

+ cosi>o[sin 2 P(D fl + cD s ) + cos 2 P(D,., + eDu) 
+ 2 sin P cos P(D M + eD w )] 

and 

C 3 = / sin co cos co sin ij(l + c cos v) dv 

= sin ia,[(cos a P - sin 2 P)(D 13 + eDw) 9) 

+ sin P cos P (A, - D 9 + eD 7 - eDu)] 
+ cos ijo[(cos 2 P - sin 2 P) (D, 4 + cD,«) 
+ sin P cos P (D 6 - Dio + eD a - eDw)] . 

Since 

E = / U. cos 7,(1 + e cos w) rfw ( 94 ) 

Jo 

we expand as with C, to yield 

E - UfE, + M^i + (Ut + JiW (123) 
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so that 

Ei — I cos 2 o) cos ??(1 + e cos v) dv (124) 

- cos ^[cos 2 P(D b + eD,) + sin 2 P(D 9 + eD u ) 

- 2 sin P cos P(D n + eDi 8 )] 

- sin T/olcos 2 P(D 6 + e/) 8 ) + sin 2 P{D W + C Z) 12 ) 

- 2 sin P cos 7 J (Z) 14 + e D 16 )] 
and 



A' 2 = / sin 2 w cos 77(1 + e cos y) 



= cos ^[sin 2 P(D b + eD 7 ) + cos 2 P(Z> 9 + eD u ) 

+ 2 sin P cos P(D 13 + eA 6 )] (127) 

- sin jjo[sin 2 P(Z) 6 + eD 8 ) + cos 2 P(Ao + e/J 12 ) 
+ 2 sin P cos P{D U + el> 16 )] 
and 



fir 

E 



3 = / sin oj cos cu cos »(1 -4- e cos y) rfy 

= cos tjoKcos 2 P - sin 2 P){D U + eZ^) 

+ sin P cos P(D 5 - A + eD 7 - eD n )] 
- sin 7jo[(cos 2 P - sin 2 P) (D u + e£) 16 ) 
+ sin P cos P(Z) 6 - Ao + eD 8 - eD l2 )\. 



APPENDIX D 



(128) 



(129) 



Case II Expansion 

We will now evaluate (91), (92), (93) and (94) once again in a manner 
similar to Appendix C, but this time using the approximation 

M = v - (A/6) siny. (102) 



(143) 
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The same notation will be used for the integrals 

A = I sin 77(1 + c cos v) dv 
J a 

= sin tjo / jcos bv cos (X sin v) + sin bv sin (X sin v) j 

• (1 + e cos v) dv 
+ cos 770 / | sin bv cos (X sin v) — cos bv sin (X sin v) \ 

• (1 + e cos v) dv . 

I'sing the approximations of (107) and (108), we have 

,1 = sin t,„[D, + clh + X(D 17 + eDu) - (X 2 /2)(D 9 + eD n ) 

- (X 8 /6)(Z)i 8 4-eD w )] 

(144) 
4- cos 7jo[i>o + cD A - X(D M + eD 18 ) - (X-/2) 

(D 10 + eDu) + (X 3 /6)(D 21 + eDu)]. 
In like manner, 

# = / cos 77(1 + e cos y) dv (145) 

= cos 7j [D, + eD 8 + X(D 17 + eD u ) - (X 2 /2)(D 9 + eD u ) 

- (xV6)(Dia+ eD u )] 

(146) 

- sin t/ [D 2 + eD 4 - X(D 20 + eD w ) - (X72)(D 10 

+ aDn) + (X 3 /o)(D 21 + eD»)l. 
Similarly, 

Ci = / cos" w sin 77(1 + e cos v) dv (147) 

C, - sin t,u[cos 2 P\(D & + eD 7 ) + X(D 16 + eD 26 ) - (X 2 /2)(D 23 
+ eDa) - (X 3 /6)(D 3 , + eD M )j 4- sin 2 P{(D 9 + eD u ) 
+ X(D M + eD 19 ) - (X 2 /2)(D 30 + eD M ) - (X 3 /6) 

• (Dss + cD 4 i)} - 2 sin P cos P{ (D 13 + eD 1B ) 4- X(D 12 
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+ eD 2i ) - (X 2 /2)(Z) 22 + eD 32 ) - (\ 3 /6)(D 37 + e/) M )}] 

(148) 
+ cos 7/o|_cos 2 P{ (Z) 6 + eD s ) - \(D ib + cD 28 ) - (X 2 /2) 

•(D 24 + eD 27 ) + (\76)(£>32 + eD 3i )} + sin 2 P{(D l0 
+ eDu) - \(D 2l + eD 22 ) - (X 2 /2)(Z> 29 + cD 31 ) 
+ (X a /6)(D M + eD 42 )} - 2 sin P cos P{(D U + eD M ) 

- \(D n + eD M ) - (X 2 /2)(Z> 19 + eD 81 ) + (XY6)(2> M 

+ flZ>«)}]> 

To simplify the writing of the mean magnetic torque equation, let 
F, = (D B + r/V) + X(Z) I6 + eDn) - (X 2 /2)(7) 23 + e£> 26 ) 

- (X 8 /6)(Dsi + eD 3s ) 

ft = (D 9 + eDu) + X(Dis + eDn) - (\*/2)(D 30 + eDas) 

- (tf/Wu + eDu) 

F 3 = (D u + cD a ) + X(D, 2 + eD 2i ) - (X 2 /2)(D 22 + eD 82 ) 

(124) 

- (x76)(D» + eZ>») 

F A = (A, + el> 8 ) - X(D 15 + cD ffl ) - (X 2 /2)(Z) 24 + eD 27 ) 

+ (X ;, /G)(Z>3 2 + ^34) 
F 8 = (D 10 + eD^) - X(D 21 + e7) 22 ) - (X 2 /2)(D 29 + eD 31 ) 

+ (X 3 /0)(D36 + eD 42 ) 
F« = (Dm + eDi 6 ) - X(D U + eD 23 ) - (X 2 /2)(D 19 + eD 31 ) 

+ (\ 3 /G)(D 3fi + eD i0 ). 

By the similarity of C 2 and C 3 integrals to C\ and of E\ , E 2 , E 3 to Ci , 
C 2 , C 3 we immediately write and summarize the following: 

Ci = sin 7j [Fi cos 2 P + F 2 sin 2 P - 2F 3 sin P cos P] 

+ cos y [F 4 cos 2 P + P 5 sin 2 P - 2F 6 sin P cos P] 

C 2 = sin T)„[Fi sin 2 P + F 2 cos 2 P + 2P 3 sin P cos P] 

+ cos tjo[P4 sin 2 P + F 5 cos 2 P + 2F 6 sin P cos P] 

C 3 = sin tjo[(Fi - F 2 ) sin P cos P + F 3 (cos 2 P - sin 2 P)\ 

+ cos t? [(F 4 - P B ) sin P cos P + F 6 (cos 2 P - sin 2 P)] 



(128) 
(129) 
(130) 
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Ei = cos tjoIFi cos 2 P + F 2 sin 2 P - 2F 3 sin P cos P] 

- sin 77o[F 4 cos 2 P 4- F 5 sin 2 P - 2F 6 sin P cos P] 

F, = cos ij [Pi sin 2 P + F 2 cos 2 P + 2F 3 sin P cos P] 

o •> ,„ . ,-, „, (132) 



sin va[Fi sin 2 P + F 5 cos 2 P 4- 2F 6 sin P cos P] 
?o[(Fi - Ft) sin P cos P + F 3 (cos 2 P - sin 2 P)] 
sin rjo[(P4 - F B ) sin P cos P + F 6 (cos 2 P - sin 2 P)]. 



— oui ';(IL« 4 o'« * i*o ^^-^ * i — ■ o - — ■ — - j 

F 3 = cos tjo[(Fi - Ft) sin P cos P + F ;! (cos 2 P - sin 2 P)] 

V loo J 



APPENDIX E 



The D Integrals'* 



f 1 

Di = / cos (to) dv = T sin (2irb) 
Jo o 

D 2 = [ sin (bv) dv = J [1 - cos (2ttd)] 

Z) 3 = / cos (bv) cos 



v dv = — 



b sin (2tt6) 



1 - & 2 
D 



= j sin (to) cos v dv = — - — ^ [1 — cos (2wb)] 

Jo J- "" 

D 5 = I cos (to) cos 2 v dv - " _ , 2 , sin (2tto) 



D 6 = J sin (bv) cos 2 w dv = w _ . [1 - cos (2ir6)] 
D 7 = /" cos (6.) cos 3 B dv = - — ^j^L_. sin (2*6) 

D 8 = £' sin (to) cos 3 V <fo = - d^y^'ly) t 1 " cos < 2 * 6 > ] 

2 sin (2tt6) 



D 9 = I cos (to) sin 2 v dv = 
Jo 



6(4 - 6 2 ) 



/• 2 * 2 

Dio = / sin (to) sin 2 a dv = ,. - , 2 , [1 — cos (2irb)] 
Jo 0(4 — ) 

f- T ■ i n ^ i 2o sin ( 27rD ) 

l^ii = / sm v cos v cos {bv) dv — — _ , 2 wq _ fo2-\ 
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D u = J sin 2 v cos v sin (bv) dv = j—————— [1 - cos (2tt6)] 

n f' W ■ n. \ j 1 "" cos (2x6) 

L>i3 = / sin v cos y cos {bv) dv = ^ — - 

Jo 4 — o z 



'0 
.2* 



,. f' r . . ,, v , sin (2tt6) 

L>i4 = / sin v cos v sin (6y) ay = — — — — 

Jo 4 — 6 2 



D» = J sin v cos 2 v cos (6w) dv = _' fo2 wi _ ^S t 1 ~~ cos ( 2r k)] 



'o x ' (9 - 6 2 )(1 - 6 2 ) 

•2t /o ,.2>, 

'Hi _ 



,, [" • 2 • /. \ 7 (3 — 6 ) sin (2x6) 
/>i 6 = / sin w cos v sin (6y) dy = - — ; 7 

Jo (9 — er)(l — 6 J ) 

,, f " . ,, <. , sin (2x6) 

L>n = / sin y sin (bv) dv = — — — — 

Jo 1 — 6 2 

r, f * • 3 . /i \ , 6 sin (2^6) 

D U = y t sm . sin (to) * = - (9 _ y)(1 _ y) 

f r . 3 • /l \ i 6 sin (2x6) 

= / sin y cos y sm (6y) ay = — 

Jo (lfi - 



/; 



o (16 - 6 2 )(4 - 6 2 ) 

n f ' ■ (L \ j 1 — cos (2<r6) 

D w = I sin y cos (yy) ay = — — — 

Jo \ — b l 

n f * • 3 /; \ j 6[1 — cos (2x6)] 

D a = I sm . co S (to) <fo = (9 _ 6!)(1 _ W) 

n /" * • 3 /l N j 6[1 — COS (2x6)] 

D, = ^ sm . cos . cos (b„) dv = (16 - y)(4 _ ^ 

n [ 2T . 2 2 /i \ j 2 sin (2x6) 

l>23 = / sin y cos v cos (by) ay = — —f 

Jo 6(16 — 6 2 ) 

n /* * • 2 2 /i \ j 2[1 — cos (2x6)] 

/A, = / sin y cos y sin (6y) dv = — — — — - — 

Jo 6(16 — b l ) 

n f- T .2 3 /, v . 26(13 - 6 2 ) sin (2tt6) 

D. = / sm y cos v cos (bv) dv = - _ _ _ 



o (25 - 6 2 )(9 - 6 2 )(1 - 6 2 ) 



D 



f 2 " 
= I sin v cos 3 v sin (6y) dv = — 

Jo 



3 • /, x , (10 - 6 2 ) sin (2x6) 

i,m(t,) *"" (iR)B 

•2ir , / 1r> ,2> 



/o (16 - 6 2 )(4 - 6 2 ) 

n Z" 2 ' . 2 3 „ ■ /, v , 26(13 - 6 2 )[1 - cos (2tt6)] 

D 27 = ] sin v cos v sm (bv) dv = - _ _ _ 



n [ • 3 ,, w (10 - 6 2 )[1 - cos (2tt6)] 

D 28 = sin y cos v cos (6y) rfy = 7 — =—71 r^ 

.'0 (16 — 6 2 )(4 — b 2 ) 
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n f' 1 ' ■ 4 ■ n \ i 24[1 " cos (2x6)] 

I>29 = / sin v sin (bv) dv = —\ — — — 

Jo 6(16 - 6 2 )(4 - 6 2 ) 

n f 2 ' ■ 4 /, v . 24 sin (2x6) 

L)-io = I sin v cos (bv) dv — 
Jo 

D 



o 6(16 - 6 2 )(4 - 6 2 ) 

f -a 2 /, x , 6(5 — 6 ) sin (2x6) 

ai = / sin v cos /> sm (6/ 1 ) dv — — : -^ 

Jo (25 - 6 2 )(9 - 6 2 )(1 - 6 2 ) 

rx f * • 3 2 /i \ j 6(5 — 6 2 )[1 — cos (2x6)] 

D„ = j o sm . cos . cos (to) * = (25 _ y)(9 _ y)(1 _ y) 

n [ 2r ■ 3 3 • /, v , 6 sin (27r6) 

O js = I sm , cos „m (to) dv = ~ (36 _ y)(4 _ y) 

n f 2 ' -s 3 /, v , 6[1 — cos (2?r6)] 

D 34 = J q sin , cos v cos (to) A = (36 _ y)(4 _ y) 



n /* -8 • /, x , 120 sin (2x6) 

D 35 = I sin , sin (bv) dv = - (2fi _ y)(p _ __ 



/• 2T • B „ w 12 

36 — / sin y cos (6v) tft; = — — — 

Jo (25 — 



Air, - 



120[1 - cos (2tt6)] 



(25 - 6 2 )(9 - 6 2 )(1 - 6 2 ) 

n f U ■ 4 • ,u\ j 246[1 - cos (2x6) J 

I> d7 = / sin v cos v sin (bv) dv = — — — -^— -— — — — 

Jo (25 - 6 2 )(9 - 6 2 )(1 - 6 2 ) 

n /" * • 4 /; \ » 246 sin (2x6) 

L>38 = / sin v cos i/co3 (6r) rfy = — — — — — - 

Jo (25 - 6 2 )(9 - 6 2 )(1 - 6 2 ) 

4 .->,,,, 24(6 - 6 2 )[1 - cos (2x6)] 



Aw = / sin 4 v cos"' v sin (6c) dv — 
J ii 



6(36 - 6 2 )(16 - 6 2 )(4 - 6 2 ) 

n f 2 ' -4 2 n \ i 24(6 — 6 2 ) sin (2x6) 

D i0 = sin v cos v cos (bv) dv = — — l2W1f . rrvri i^ 

J« 6(36 — 6 2 )(16 — 6 2 )(4 — 6 2 ) 

120 sin (2x6) 



D41 = / sin"' ;■ cos v sin (bv) dv = — 

Jo (36 - 6 2 )(16 - 6 2 )(4 - 6 2 ) 

n C .6 ,, x , 120[1 - cos (2x6)] 

Dv> = / sin v cos v cos (6r) dv = — — r^—, 

J., (36 - 6 2 (16 - 6 2 ) 4 - 



(36 - 6 2 )(16 - 6 2 )(4 - b°-) 
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